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ABSTRACT 


A computer model simulating the seasonal variations of 
mixed layer nutrient concentrations, phytoplankton biomass 
carbon, and herbivorous zooplankton biomass carbon was devel- 
oped, The simulation was generated using an annual cycle of 
four environmental parameters: (1) incident solar radiation, 
(2) upwelling velocity, (3) mixed layer depth, and (4) mixed 
layer temperature. Simulation results were compared with 
nutrient and zooplankton biomass data collected on a series 
of seven cruises made in central Monterey Bay from February 
through December, 1974. Both observed and simulated Zoo- 
plankton stocks were characterized by two distinct maxima. 
The initial peak (1.05 gC /m”) occurred in late July and was 
followed by a decline in populations through the month of 
August. During the fall and early winter, zooplankton biomass 
increased rapidly to an overall maximum of 1.85 eC /m?. 
Individual environmental parameters were tested to ascertain 
their importance in controlling simulation results. Phyto- 
plankton stocks were influenced principally by changes in 
incident radiation, whereas temperature variations produced 
the most significant fluctuations in zooplankton biomass. 
Simulation responses suggest that upwelling, in addition to 
providing nutrients for primary production, enhances zZoo- 
plankton productivity by bringing colder deep water to the 


surface, thereby reducing zooplankton respiration requirements. 
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I, INTRODUCTION 


A. PURPOSE 

There presently exists a large body of detailed informa- 
tion about the biological, physical, and chemical processes 
which occur in the marine environment. These processes act 
in concert to control the plankton populations of the sea. 
The analysis and interpretation of these environmental and 
biological factors has been significantly enhanced by the 
development of computer simulation techniques. The aim of 
computer modeling is to reproduce the dynamic response of a 
biological ecosystem to the complex mix of interacting 
variables which act on and within that system. 

The purpose of this research was threefold: 


(1) to formulate a model of a planktonic ecosystem in a 
region of upwelling, 


(2) to develop a computer Simulation of the model which 
generates seasonal changes in mixed layer nutrient 
concentrations, phytoplankton biomass carbon, and 
zooplankton biomass carbon, and 

(3) to test the accuracy of the simulation results by 
comparison with observed data. 

B. CONCEPTUAL MODEL 

The first step in the modeling process consists of 
reducing the natural system being simulated to a simplified 
conceptual model. From this conceptual model, processes and 


rates are formulated into mathematical equations. The model 


is completed by translating the equations into the proper 
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format. The conceptual model considered here describes the 
relationships between factors controlling biological 
productivity (see Fig. 1). 

i. Prodwetion 

The synthesis of organic compounds from the inorganic 

constituents of the sea by the activity of living organisms 
is termed "production" (Tait and DeSanto, 1972). This 
Synthesis is carried out almost exclusively by the photo- 
synthetic activity of marine plants. The raw materials for 
photosynthesis are water, carbon dioxide and various nutrients, 
mainly inorganic ions of nitrogen (NO, NO,, NH) and 
phosphorus (PO, HPO,, HyPO;). 


plants utilize light energy (nhy) to synthesize higher energy 


The chlorophyll-containing 


organic compounds from the lower energy inorganics as follows: 


6nhy + 6CO, + 6H,O0 + 6CH,O + 60, (1) 


This synthesis is referred to as "gross primary production." 

A certain percentage of this synthesized material is 
broken down by the plant cells to provide energy to carry on 
physiological processes. This loss is termed respiration. 
The remainder of the synthesized material iS available for 
cell growth and replication and is referred to as "net 
primary production." It is the net primary production of 
phytoplankton that is available for consumption by the 
herbivorous zooplankton, 

Grazing by herbivores on the phytoplankton leads to 


"Secondary production."' Herbivores in turn are a food source 
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for first level carnivores ("tertiary production"), who them- 
selves become prey for other carnivores, and so on. This 
type of organization of a biological community is referred 
to as a "trophic structure" or “food chain," At each level 
losses of biomass may occur from respiration, excretion, 
predation, or natural mortality, In reality the relation- 
ships are more complicated than this simple "food chain" 
concept indicates. Since one organism may take food from 
more than one trophic level, a more correct structure would 
define an intricate “food web.™ 

In spite of its over-cSimplification, the "food chain" 
model depicts certain important characteristics of the marine 
planktonic ecosystem. First, the path of energy in the form 
of organic material is part of a closed cycle in which 
mucvrrentSware assimilated, transferred, and regenerated; and, ~ 
Since water and carbon dioxide are abundant in the sea, the 
factors which potentially limit photosynthesis are the 
availability of nutrients and light energy to the phytoplankton 
(Ryther, 1963), and the temperature of the medium (Parsons 
and Takahashi, 1973). | 

2. Light 

The light energy available to the phytoplankton is 
dependent on two parameters, the amount of radiation incident 
on the sea surface (usually expressed as g-cal/em“/min) and 
the transparency of the water (expressed as the extinction 
coefficient, k). Ryther (1956) observed that photosynthesis 


increases linearly with light intensity up to a saturation 
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point. Beyond this level greater illumination produces no 
further increase in photosynthesis and eventually results in 
an inhibition of the photosynthetic rate. 

Light penetrating into the water is absorbed and 
scattered selectively. The intensity is attenuated exponen- 
tially with depth and can be expressed in terms of the 


extinction coefficient (k): 


as -KzZ 
I. ~ Ie (2) 
I, = radient energy at the surface 
I, = radient energy at depth z 
k = light extinction per meter depth. 


The transparency of the water is dependent on the 
concentration of phytoplankton and the presence of particulate 
or dissolved substances which absorb light. If we disregard 
selective absorption of light of different wavelengths and 
consider the total visible spectrum, we may say that for 
clear ocean water, k equals 0,04 - 0.05 (Ryther, 1963). Riley 
(1956) determined empirically the relationship between 
phytoplankton density and transparency. He assumed that 
absorption by material other than water and plants was 


negligible. His relationship is: 


k = 0.04 + 0.0088 C1 + 0,054 12/9 (3) 
where Cl is chlorophyll concentration (yug/1) and 0.04 is 
the extinction coefficient for pure sea water. It is apparent 


from this formula that one mechanism for limiting phtoplankton 
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y 


productivity is self-shading. Higher concentrations of algae 
produce a higher extinction coefficient and, from equation (2) 
reduce the light intensity in the water column. The reduced 
intensity lowers the rate of photosynthesis. 

Measurements of the photosynthetic rate under 
conditions of bright daylight indicate that radiation intensity 
near the surface is at or above the saturation point for most 
phytoplankton. Under these conditions maximum productivity 
occurs at a depth, usually between 5 - 20 m. Below the depth 
of maximum productivity the photosynthetic rate decreases as 
the light intensity is attenuated. A point is reached where 
the rate of photosynthesis is just sufficient to balance the 
rate of breakdown of organic materials by respiration. This 
point is referred to as the "t'compensation depth" (Tait and 
Besanto, 1972). 

A second level below the compensation depth is the 
"critical depth."' The concept of the "critical depth" was 
first suggested by Grann and Braarud (1935) and developed into 
a mathematical model by Sverdrup (1953). The critical depth 
is the level at which the total respiration of the phyto- 
plankton in the water column just equals the total photo- 
synthesis. The distance between the compensation depth and 
the critical depth depends on the proportion of the 
phytoplankton stock above and below the compensation depth. 

If the standing stock of phytoplankton is to increase, 
its total phtosynthesis must be greater than its total 


respiration, This is possible as long as the depth to which 
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the plankton is mixed does not exceed the critical depth, 
When the depth of the mixed layer extends below the critical 
depth (assuming uniform plankton mixing), total losses exceed 
total gains and the stock declines. 
3. Nutrients 

Another parameter which can critically limit the 
growth of phytoplankton is the concentration of inorganic 
nutrients in the upper, "euphotic" layers of the ocean. Many 
elements such as manganese, silicon, copper, molybdenum, and 
cobalt are essential for the growth of particular species. 
However, ions of nitrogen and phosphorous may be present in 
sufficiently low concentrations that their availability exerts 
a dominant control over production (Tait and DeSanto, 1972). 

The absorption of nutrients in the surface layer by 
phytoplankton can reduce the concentration of these ions to 
levels where further uptake is inhibited. The kinetics of 
this inhibition was investigated by Ketchum (1939). He 


demonstrated that the uptake of nitrate by the diatom 


Nitzschia closterium was concentration dependent over the 


approximate range of 1 - 7 ug-atom N/1. He also showed that 
phosphate followed a similar inhibition. Dugdale (1967) 
described the rate of nutrient uptaking using Michaelis- 


Menten enzyme kinetics in which: 


mm N 
i ene Ko + N ee) 


v = rate of nutrient uptake 
V = maximum rate of nutrient uptake 


Max 
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K. = half saturation constant, i1,e. the nutrient 
concentration at which v = V Vie 
max 
N = nutrient concentration 


Studies conducted by Eppley and Thomas (1969) on two species 
of diatoms indicate that the half-saturation constants for 

growth and nutrient uptake are very Similar. This means that 
the expression given for velocity of nutrient uptake can also 


be used to describe the growth rate: 


- N 
H = max KO +N sl 


ll 


Lu growth rate 


maximum growth rate 
By plotting N/u versus N, a straight line is obtained with 
the intercept on the abscissa at “K.. 
a. Temperature 
Ichimura (cited in Parsons and Takahashi, 1973) 
studied the relationships between phosphate concentration 
and the rates of organic synthesis in natural populations of 
phytoplankton in Tokyo Bay. His results indicate a definite 
Michaelis-Menten type response which is highly temperature 
dependent (see Fig. 2). 
b. Excretion 
A certain fraction of the nutrients absorbed by 
phytoplankton are returned to the medium both in the euphotic 
zone and at depth by zooplankton. According to Harvey (1957), 
nearly all phytoplankton is grazed by herbivorous zooplankton. 
A portion of this grazed algae can be voided in an undigested 
or semidigested state. It is unclear, however, how much 


phytoplankton nutrient elimination by zooplankton is returned 
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in an available form in the euphotic zone and how much sinks 
to greater depths before being remineralized. Some authors 
(Gardner, 1937; Riley, 1951; and Harris, 1959) consider 
zooplankton to be a Significant source of useable nutrient 
am the euphotic zone, while others believe this source to be 
negligible (Steele, 1959; Rigler, 1971). Most recent work 
seems to favor the higher estimates. A distribution of 
dietary phosphorus in Calanus during active feeding (April) 
has been given by Butler et al. (1970) as 17.2% retained for 
growth, 23.0% voided as fecal pellets and 59.8% excreted as 
soluble phosphorus (organic and inorganic). These estimates 
indicate that over 80% of the phosphorus ingested by Calanus 
gets back into the environment and nearly 60% is soluble 
phosphorus. The ability of marine algae to utilize soluble 
organic phosphorus was demonstrated by Kuenzler (1965). 

The importance of zooplankton as a source of 
phytoplankton nutrients was verified in two studies cited by 
Parsons and Takahashi (1973). Cushing (P. 139) "showed that 
during ten weeks of the spring phytoplankton/zooplankton bloom 
in the North Sea, inorganic phosphate did not decrease below 
0.6 wg-atom/1 due to zooplankton excretion. However, Antia 
et al. (1963) observed a decrease in phosphate to 0.1 ug-atom/1 
after two weeks of phytoplankton growth in the absence of 
zooplankton," suggesting the important role of zooplankton as 
a source of phytoplankton nutrients. 

The nitrogen and phosphorous cycles in the sea 


follow Mmsically the same course, However, the remineralization 
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of nitrogen occurs at a much slower rate than phosphorus. 

In the experiment conducted by Antia et al. (1963) half the 
phytoplankton phosphorus was remineralized in 14 days, but 
after 75 days no remineralization of nitrogen was observed. 

A second major difference between the cycles is 
that all three inorganic nitrogen compounds, ammonium (NH), 
nitrite (NO,), and nitrate (NO) can be absorbed by phyto- 
plankton. Absorption is not indiscriminate, however. There 
is a marked preferential uptake of ammonium (Dugdale and 
Goering, 1967). The nitrogen cycle then consists of three 
seperate cycles, one for each of the three nitrogen containing 
ions involved in photosynthesis, and possibly a fourth involv- 
ing urea (Dugdale, personal communication). 

Although a fraction of the nutrients absorbed by 
the phytoplankton during photosynthesis are recycled within 
the euphotic zone, the remainder are regenerated in deeper 
water. The result is a net depletion of nutrients from 
surface waters and an accumulation of nutrients at depth. 

The continuation of plant growth in the sea, therefore, is 
largely Res ater upon the various processes which return 
nutrient-rich water from below the euphotic zone. 

c. Mixing and Upwelling 

Turbulent (eddy) diffusion is a primary mechanism 
of injecting nutrients into surface waters of many ocean areas. 
Vertiché eddies may be generated by a variety of dynamic 
forces including thermal convection, internal waves, friction 


at current boundaries, and wind-mixing. Turbulence may have 
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a variety of effects on biological production, It may 
enhance production by bringing nutrients to the surface, or 
it may reduce production by carrying phytoplankton down 
below the euphotic zone. 

A second major phenomenon operating in more 
localized regions is coastal upwelling. For example, along 
western continental boundaries in the Northern Hemisphere a 
northerly wind stress blowing parallel to the coast causes 
an offshore transport of water with a compensating replacement 
of upwelled waters from 200-300 m depth (Sverdrup et al., 
1942), 

An index of the intensity of upwelling has been 
estimated using wind stress data (Wooster and Reid, 1963; 
Bakun, 1973). This “Ekman transport" of water away from the 


coast is calculated according to the following equation: 


M = woe (6) 
M = "hkman transport" normal to the coast 
i = Wind stress parallel to the coast 
f = Coriolis parameter 


Using the magnitude of offshore transport it is 
possible to compute the velocity of the associated verticle 
motion. Wooster and Reid (1963) assumed a steady-state 
offshore "Ekman transport" of 10 kg/cm/sec with a coastal 
upwelling band 50 km wide and estimated the compensating 
vertical motion to be about 50 m/month, a resut compatible 


with previous estimates which ranged from 10-20 m/month to 
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80 m/month (Smith, 1968). Although the depth from which these 
waters are upwelled is not believed to exceed a few hundred 
meters, the input of nutrients is sufficient to maintain 
the regions of upwelling at a high level of productivity. 
Horizontal advection of surface waters could have 
either a positive or negative effect on nutrient concentration. 
If “upstream" or “offshore waters have been depleted below 
teeat Values, advection into the area under consideration 
would reduce the concentration and conversely, higher values 
would raise local concentrations. This process is discussed 
Paver . 
4. Grazing 

Another major factor controlling the productivity of 
phytoplankton is grazing by herbivorous zooplankton. Accord- 
ing to a number of authors (Steemann-Nielsen, 1958; Cushing, 
1959; Steele, 1974) grazing by zooplankton is the primary 
factor affecting the dynamics of plankton populations. 

Mullin (1963) identified two methods of describing 
Poeoptlankton grazing activity, by the "filtering rate" and 
the "food intake rate." The "filtering rate" is the rate at 
which water is swept clear of food organisms. The "rate of 
food intake" is a measure of the biomass of food organisms 
ingested per unit time and depends on the filtering rate and 
the concentration and size of the food organisms. Mullin 
(1963) demonstrated that, for four species of copepods, the 
volume swept clear is inversely proportional to food concen- 


tration. Steele (1974) showed that such an inverse relation 
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corresponds to a hyperbolic rate of intake of biomass. As 
the food concentration increases, the rate of biomass 
ingestion approaches a maximum (see Fig. 3), McAllister (1970) 
observed intake rates which approximate the type given by 
Steele, Wash and Bass (1971) expressed this saturation feed- 
ing behavior in terms of a Michaelis-Menten equation. This 
hyperbolic response suggests two important characteristics of 
herbivore grazing. First, a minimum phytoplankton concentra- 
tion exists at which grazing ceases; and second, as algal 
concentration increases, herbivores asymptotically approach 
a Maximum ingestion rate. 

In summary, it iS apparent that the productivity of 
a region is dependent on the availability of light and 
nutrients to the phytoplankton, the temperature of the medium, 
and the grazing and excretion of biomass and nutrients by 
the herbivores. The rates at which these processes proceed 
are not simple functions but are dependent on a complex mix 
of physical, chemical, and biological parameters. In order 
to model the dynamics of this regional plankton community, 
therefore, it is essential to utilize a method of computer 
Simulation. Simulation provides a technique for integrating 
the many variables affecting productivity into a system of 


interacting components, 
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II, BACKGROUND THEORY 


A. EARLY WORK 

The basic equations of biological growth kinetics were 
developed in the early 1920's by A. J, Lotka (1925). 
Volterra (1926) extended Lotkats equations to describe predator- 
prey interactions. The Lotka-Volterra equations (7) consist 
of a pair of differential equations in which the growth of 
each population is limited by two factors, self inhibition 


(ry, ro) and competition (a, b) (Odum, 1971): 


K, - N, «aN 


GINe se 1 i 2 
de> <b iam . bey 
aN, a. Ko - No _ DN, 
dt Zee, Ky 
Ny, No = quantity of species one and two 
a, b = competition coefficient indicating the negative 
effect one species has on the other 
ry, To = the intrinsic growth rates of each species 
Ky» Ko = carrying capacity, the largest population of 


each species the environment can support. 

Fleming (1939) used a modified form of the LotKa- 
Volterra equations to calculate the seasonal change in 
phytoplankton population: 


dP 


ape = Gs eta) : ey 
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Pp = Integral phytoplankton population under a unit 
of water surface 


a = specific growth rate of plankton 

b = initial grazing rate | 
(b+ct) = grazing rate as a function of time 
Cc = rate of grazing rate change 

t = time. 


His model emphasizes grazing as the factor limiting phyto- 
plankton populations, Patten (1968) pointed out that, 
although the model may be unrealistic in specifying a constant 
growth rate and a linearly increasing grazing rate, the model 
"worked" because it accurately corresponded to observations 


made by Harvey et al. (1935) in the English Channel. 


B. PHYTOPLANKTON MODEL 

Riley (1946) made a Significant improvement upon the lack 
of realism in Fleming's model. He defined the change in 
phytoplankton as: 


dP 


apo PCP. =) Ree G) (9) 


P = integral phytoplankton under a Square unit of 
waberesuriace . 
P, = photosynthetic rate 
R = respiration rate 
G = grazing rate. 
Instead of using a constant photosynthetic rate, Riley 
fewenea Po as a function of light intensity, nutrient concen- 


h 


tration and verticle turbulence. 
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aKZ 


Ph = (pI. /kz,) (1 -e mece -~ N) (i ~~ V) (10) 
P = fraction of light used in photosynthesis 
a = surface light intensity 
k = extinction coefficient 
Z, = depth of euphotic zone 
N = nutrient depletion factor 
V = factor expressing the effect of turbulence 


His expression for the respiration rate was defined in terms 


of temperature, 
afk 


R = Roe CE) 
R = respiration rate 
Ry = respiration rate at O°C 
r = rate of respiration rate change 
T = temperature. 


The grazing rate was assumed to be proportional to the inte- 


gral quantity of zooplankton, 


G= g¢Z CizZ 


G = grazing rate 


g = proportionality constant between grazing rate and 
zooplankton 


Z= integral zooplankton under a square unit of water 
surface. 
C, PHYTOPLANKTON, NUTRIENT MODEL 
Steele (1958) introduced a model based largely on the 
equations developed by Riley. His contribution is significant 


in that it is one of the first "systems" approaches to plankton 
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dynamics. Where Riley treated nutrient concentration as an 


independent parameter, Steele defined nutrients and phyto- 


plankton as two inter-«dependent variables, 


To develop a nutrient equation, Steele made certain 


simplifying assumptions about the environment. The ocean 


volume was divided into two thermally-stratified layers with 


homogeneous mixing in each layer. Mixing across the boundary 


was controlled by an exchange coefficient (m). 


The system of equations devised by Steele define the rate 


of change of nutrients and phytoplankton in a surface layer 


40m deep: 


il 


= = (Py - rene fZ)P- + - m(P - a, (13) 
dp _ 
mae -C(p), -~ r)P -~- m(p - P,) (14) 
mixing coefficient governing exchange between the 
upper and lower layers. 
DovTommruccon concentravion in the upper layer. 
phytoplankton concentration in the lower layer. 
photosynthetic rate 
phytoplankton respiration rate 
grazing rate, a function of zooplankton populations 
zooplankton population 
ratio of sea water viscosity at O°C to that at 
ambient temperature 
factor for converting carbon to phosphate 
nutrient concentration in upper layer 
HUtrIent= concentration im lower layer. 
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D, REGIONAL STEADY STATE MODEL 

Riley (1965) added to Steele's (1958) model a term 
expressing the regeneration of nutrients voided by zoo- 
plankton. This model was tested against observations made 
during quasi-steady state summer conditions found in various 
regions of the Atlantic and Pacific Oceans. Basic environ- 
mental variables, including phosphate concentration in the 
deep water mass, the rate of mixing between the two layers, 
and water temperature, were varied according to the charac- 
teristics of each area. For each array of environmental 
factors, steady-state concentrations of phytoplankton, 
zooplankton, and phosphate were calculated for the upper 
layer. Riley concluded that the model "provided a reason- 


able fit for most of the observed variations," 


E. SHORT TERM SPATIAL MODEL 

Walsh and Bass (1971) developed a spatial model of an 
upwelling ecosystem which was used for short term analysis. 
Their system of equations related the upwelling of nutrients 
to the nitrogen biomass of phytoplankton, zooplankton, and 
fish. A Michaelis-Menten type expression was incorporated 
into the model to define both the nutrient uptake rate and 
the grazing rate. 

To conclude, the early models of Riley (1946) and Steele 
(1958) established the fundamental principles used in 
Simulating plankton dynamics. Riley (1965) demonstrated 


that regional variations in phytoplankton could be 
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realistically predicted for a quasi-steady state summer 
period. Using local environmental conditions as input 
parameters, his model accurately predicted the plankton 
stocks observed in different geographic areas, Walsh and 
Bass (1971) developed a spatial model to analyze a specific 
region of upwelling. Since their model applied to short 
term (two week) simulations, environmental parameters 


were held constant. 
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Ill, THE MODEL 


The simulation model proposed in this thesis defines a 
one dimensional, nonlinear, planktonic ecosystem in a 
centralized region of coastal upwelling (Monterey Bay). 
The simulation output consists of a theoretical annual cycle 
of mixed layer nutrient concentration, phytoplankton and 
zooplankton biomass. Although the model has been developed 
for a localized area it may be considered applicable to a 
broader geographic region, Bolin and Abbott (1962) determined 
that "results obtained (for the Monterey Bay area) should 
reflect in a broad way the oceanographic conditions and 
events occurring over an extensive segment of the eastern 
North Pacific."* As in Riley's (1965) model, environmental 
conditions control the output of the system. From observed 
Semanezes in the physical environment, the model generates 
values for phosphate concentration and plankton standing 


stocks. 


A, COMPARTMENT MODEL 

The compartment model (Fig, 1) defines the flow of 
nutrients and organic carbon through a three component eco- 
system consisting of inorganic phosphate, phytoplankton, and 
herbivorous zooplankton, In this diagram (method from Patten, 
(1971)) the state variables of the system are denoted by 


blocks, and signal flows between the "compartments" are 
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designated by unidirectional branches, The relations between 
the state variables are expressed as a system of differential 
equations, 

The three equations in Fig, (1) describe the change in 
mnosphate concentration in the mixed layer (X11), in the bio- 
mass of phytoplankton (X2), and in the biomass of herbivorous 
zooplankton (X3). The underlying principle used in the 
derivation of these expressions is the conservation of mass, 
Therefore, the validity of the system of equations depends 
on having all significant transfers of matter accounted for. 

The movement of substance or energy into or through the 
system is termed a "flux" or "flow" and is described as a 
quantity of mass or energy per unit time, e.g., gram calories 
per square centimeter per day or grams carbon per square 
meter per day. In ecological systems the flux between 
compartments is often dependent upon the state variables 
mvolved. In such cases the flux is defined in terms of a 
"rate of transfer" or "flow rate." A flow rate is a 
quantity transferred per unit quantity of source or receiver, 
per unit time. The distinction between flux and flow rate 
can best be made in terms of units. If a flux is in grams 
Carbon per square meter per day, the corresponding rate of 
flow 1S in units of reciprocal days. 

To differentiate between fluxes and flow rates in the 
model, fluxes have been denoted by descriptive variable names 
(NUT, REGEN, GRAZ, etc.). Flow rates are designated by a 


Greek letter followed by two digits (TAU12, RHO30, LAM30, etc.). 
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“TAUY identifieswa transfer between tropic levels, "RHO" 
signifies respiration, and "LAM" represents losses due to 
mortality and carnivorous predation. The first digit follow- 
ing the Greek letter identifies the source compartment and 
the second, the receiver. A zero indicates a flow into or 


out of the system. 


B. ECOSYSTEM INPUTS AND BOUNDARY CONDITIONS 

Inputs which originate as energy or material sources 
outside the system act as forcing functions. The forcing 
functions which drive the system of equations consist of 
muirecnvironmental parameters, incident solar radiation (RADI), 
PemeoUriace Lemperature (TEMP), a coefficient of water volume 
exchange across the layer (M), and an estimate of upwelling 
velocity (CW). 

Values for incident solar radiation (RADI) were obtained 
from tables developed by Kimball (p, 103, in Sverdrup et. al., 
1942), The values desired for latitude 36° N. were estimated 
by interpolating between those tabulated for 42° N. and 30° N. 

Temperature data was obtained during a series of seven 
cruises conducted aboard the R/V ACANIA from February through 
November 1974, Estimates of the mixed layer depth (Z), used 
fomeactine boundary conditions, were developed from the 
empirical formula of Robinson and Bauer (1971), This formula, 
currently used by Fleet Numerical Weather Central (FNWC) 
defines the mixed layer depth as "The depth at which the 
temperature of the water is 2° F, (1.11° C) below the surface 


temperature." 
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Bakun (1973) developed a series of monthly indices 
describing the intensity of large scale, wind induced coastal 
upwelling at selected locations along the west coast of 
North America. The indices are based on calculations of 
offshore Ekman transport derived from monthly mean surface 
atmospheric pressure data. Values of offshore transport for 
1974 (Bakun, Personal Communication) were used to estimate 
upwelling velocity (W), A maximum velocity of one meter 
per day was assumed to correspond with the maximum index of 
upwelling. A linear interpolation was then applied to other 
monthly indices to arrive at an annual variation of upwelling 
velocity, The coefficient of exchange (M) was assigned a 
constant value of 0,12, Other sources of nutrient input 
(i.e., river runoff etc.) have been assumed to be negligible. 

As mentioned earlier, the models of Steele (1958) and 
Riley (1965) postulated a two layered ocean with homogeneous 
mixing in each layer. Phytoplankton growth was confined to 
the upper layer and the phosphate concentration was assumed ~ 
to remain constant in the lower layer, Mixing between the 
two layers supplied nutrients to the surface. The present 
model incorporates these assumptions and includes upwelling 


as an additional source of nutrient input. 


C. ECOSYSTEM EQUATIONS 

Having specified the input signals and boundary conditions, 
it is possible to define a system of equations to generate 
the response of the state variables to variations observed in 


the physical environment. 
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ft aerubolent Lquavion 


The first equation, expressing the change in phosphate 


concentration in the mixed layer, is given by: 


oe = NUT + REGEN «= UPTAK 1B) 
mee. 
The change The input The regene i The loss due 
in nutrient Ol snier hen t eration in- to photosyn- 
concentrat-(=)from upwell<-(+)put from (-) thetic uptake 
ion with ing and zooplankton by 
time mixing excretion phytoplankton. 


The expression defining the input of nutrients from the lower 


layer is! 
NUT + (CW/EKL)*(DX1~-X1)) + (C(M/Z)* (DX1~-X1) ) CEG) 

NUT = flux of nutrients into the 

surface layer (ug atom P/1i/day) 
W = upwelling velocity (m/day ) 
EKL = depth of Ekman layer (m ) 
DX1 = nutrient concentration in 

the deep layer (ug atom P/1) 
X1 = nutrient concentration in 

the mixed layer (ug atom P/1) 
M = coefficient of exchange (m/day ) 
Z, = mixed layer depth (m), 


The ratio of the coefficient of exchange to the mixed layer 
depth (M/Z) represents the fraction of the upper layer 
transferred to the lower layer in a days time, with a 
corresponding replacement from below. The ratio of the 


upwelling velocity to the depth of the Ekman layer (W/EKL) 
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represents the fraction of the surface water that is trans- 
ported offshore and replaced by upwelled water, 
The second term in the nutrient equation represents 


the input of phosphate from the nutrients voided by 


zooplankton: 
REGEN = TAU31 * X3 (17) 

REGEN = flux of nutrients from 

zooplankton (ug-atom/1/day) 
TAU3S1 = rate of input of phosphate 

into the nutrient pool 5 

by zooplankton (ug-atom/1/g C/m' /day) 
X3 = standing stock of 9 

zooplankton Ger Cc /m- ). 


The rate at which nutrients, voided by zooplankton, are made 
available for phytoplankton uptake is a function of the 


grazing rate of the zooplankton. 


TAU31 = FAVL * (C/Z) * TAU30 (18) 


TAU31 


rate of input of phosphate 
to the nutrient pool by 5 
zooplankton (ug-atom/1/g C/m /day) 


FAVL fraction of voided nu- 
trient made available 


for phytoplankton uptake (dimensionless) 


C = converts biomass to 3 
nutrients voided. (ug-atom/1/g C/m* ) 


Z = depth of the mixed 
layer (m) 


TAU3O rate at which biomass 


1s voided by -1 
zooplankton tebe) Se 
The rate at which biomass is voided is dependent upon the 


grazing rate (see equation 37). 


33 





The final term in the nutrient equation is an 


expression for the loss of phosphate to phytoplankton: 


UPTAK = C * TAUI2 * (X2/Z) (19) 


UPTAK = flux of nutrients to 
phytoplankton (ug atom P/1/day) 


e = a conversion factor 
relating organic carbon 3 
to phosphate (ug atom P/1/g C/m~ ) 
TAU12 = rate of organic pro- 
auction by a unit of 1 
phytoplankton (day) 


eZ, = standing stock of 9 
phytoplankton (g C/m™ ) 


Z, = mixed layer depth (m)., 


2. Phytoplankton Equation 

The equation for the second state variable, phyto- 
plankton, is an adaptation of Riley's (1946) model. Modi- 
fications have been made on the basis of recent findings, e.g. 
the effects of a potentially limiting nutrient on photo- 
synthesis have been expressed in terms of Michaelis-Menten 
kinetics (Dugdale 1967); McAllister (1970) demonstrated that 
zooplankton grazing could be defined in terms of saturation 
feeding; and, Walsh and Bass (1971) used a form of the 
Michaelis-Menten expression to describe both the effects of 
nutrient limitation in photosynthesis and the saturation 
response in zooplankton grazing. The approach of Walsh and 
Bass has been followed here. 

The differential equation for the change in phyto- 


plankton biomass carbon is: 
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aX2 
dt 


which in word form is: 


= PROD - RESP2 <« GRAZ (20) 


The change The proe The loss The loss to 
in phyto- duction by due to grazing by 
plankton (=) photosyn- (=) reSpira- (-) zooplankton. 
biomass thesis tion 


with time 


The expression defining primary production is a 


product of the photosynthetic rate and the phytoplankton 


biomass: 
PROD = TAU12 * X2 C21) 
PROD = photosynthetic flux to D 
phytoplankton biomass (g C/m /day) 
TAU12 = rate of organic produc-— 
tion by a unit of phyto- -1 
plankton (day ~) 
Xz = standing stock of 5 
phytoplankton (= C/m ). 


Since photosynthesis in the ocean may be limited by temperature, 
light, or nutrients, each of these parameters must appear in 


the rate equation: 


TAUI12 = P * RAD * NTLIM eZ) 


TAU12 = rate of organic production -1 
by a unit of phytoplankton (day ~) 


Pp = photosynthetic constant, 
converts light energy units 
to photosynthetic 


potential (g cal/cm?/min)7* 
RAD = mean radiant energy 

available in a water D 

column of depth (Z) (g cal/em /min) 


NTLIM 


potential photosynthetic 
rate, a function of temperature ae 
and phosphate concentration (day ~). 
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The value for the radiant intensity at depth Z (RADZ) is 


derived from the incident intensity (RADI) as follows: 


RADZ = RADI * exp(-K*Z) C23) 
RADZ = intensity at depth Z (g cal/om“/min) 
RADI = incident intensity (¢g cal/cm/min) 
K = extinction coefficient (m7) 
Z, = mixed layer depth (m ) 


To find the mean radiant energy (RAD) between the surface 
and the mixed layer depth (Z), equation eight is integrated 


mom Zero to Z and devided by Z, 


Ud 
RAD = (RADI/Z) * if exp(-K*Z )dZ 
O 


RAD (RADI /(-K*¥Z)) * (1,0 - exp(-K*Z) ) (24) 


where exp(-K*Z) = 1.0 at the surface. 


The extinction coefficient (K) was estimated using Riley's 


(1956) empirical formula: 


Kose OmOtere (O.00SS + CL) + (0.054 * CL) 2/3 (25) 
K = extinction coefficient cies 
CL = chlorophyll concentration (mg /m”) . 


Riley (1959) indicated that carbon:chlorophyll ratios may 
vary from 30 to 128. During periods of algal growth, the 
ratio approaches the lower limit, It has been assumed that 
phytoplankton densities are significant during this rapid 
growth phase. Therefore, the lower ratio of 30 was used. 
Using this ratio, the value for the chlorophyll concentration 


becomes: 
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Clee some or CX 2/Z ) . 

The function NTLIM (Fig. 2) defines the maximum 
photosynthetic rate possible for a given combination of 
temperature and phosphate concentration. The function 
approximates the results of Ichimura (p. 72, in Parson and 


Takahashi, 1971): 


NTLIM = MAXN * (X1 « XIMIN)/(KN + (X1 = XIMIN)) (27) 


NTLIM = photosynthetic rate as a 
function of nutrient con- 


centration and temperature (da 
X1 = nutrient concentration 

in the mixed layer (ug atom/1) 
XIMIN = minimum nutrient concen- 

tration below which 

phytosynthesis will not occur (ug atom/1) 
KN = half saturation constant (ig atom/ 1). 


The value for MAXN is defined as a function of temperature: 


MAXN = S * TEMP + B (28) 


MAXN = maximum rate of photo- 


synthesis at a given 1 


temperature (day ~) 
S = rate of change of MAXN re | 
with temperature (day ~/°) 
B = MAXN at O°C (“¢) 
TEMP = mixed layer temperature C2C)., 


In summary, the actual photosynthetic rate for a 
specified set of light conditions, nutrient concentration, 
and temperature is a product of the photosynthetic potential 
(P*RAD) and the maximum possible rate (NTLIM),. The second 


term in the phytoplankton equation represents the loss of 
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biomass from respiration: 


RESP2 = RHO2O * X2 


RESP2 


RHO2ZO 


X2 


(29) 


flux of biomass expended 5 

py ce li “eespiration (g C/m’ /day) 
rate of respiration per 1 

unit of phytoplankton (day) 
standing stock of phyto- 9 
plankton ee me). 


Phytoplankton respiration rate is defined as a function of 


temperature following Riley (1946): 


RHO20 = RHOZ2 * expCR*TEMP) (30) 
RHO20 = rate of respiration per 1 
unit of phytoplankton (day ~) 
RHOZ2 = rate of respiration at -1 
OG. (day ~) 
R = constant, rate of change 
in respiration rate with =a 
temperature enc: 7) 
TEMP = mixed layer temperature Ge) :. 
The final expression in the phytoplankton equation represents 
the loss of biomass to grazing by herbivorous zooplankton: 
GRAZ = TAU23 * X3 (31) 
GRAZ = flux of biomass from 


TAU23= 


X3 


phytoplankton to 


zooplankton (g C/m*/day) 
zooplankton grazing rate (day~*) 
standing stock of 9 
zooplankton CoC /m). 


A saturation grazing rate is defined by a Michaelis-Menten 


type equation (Fig. 3) following Walsh and Bass (1971): 
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TAU23 = MAXG* (X2MG-X2MIN ) /(KG+(X2MG—X2MIN) ) (32) 


TAU23 = zooplankton grazing rate (day™+) 
MAXG = maximum specific grazing rate Gen) 
X2MG = concentration of phytoplankton 3 
in a mixed layer (mg C/m- ) 
X2MIN = minimum concentration of 
phytoplankton below which 3 
Grazing wil bewWOtuOccur (mg C/m_ ) 
KG = half saturation constant (mg C/m?). 
3. Zooplankton Equation 
The zooplankton equation is! 
o** = GRAZ ~- RESP3 ~ VOID - LOSS 
which in words is! 
The change The input from The loss The loss of 
in zooplank- SlhazZ ene Oon of biomass biomass by 
ton biomass (=)phytoplankton(-) from respi- (-) excretion 


with time biomass ration 


The loss of 
biomass by 
(-—) predation and 
mortality 
The first term (GRAZ) in this equation has been defined 


(see eq. 31). The second term is a temperature dependent 


respiration term Similar to that for phytoplankton: 


RESP3 = RHO3O * X3 (34) 
RESP3 = flux of biomass during 9 
respiration (g C/m™ /day) 
RHO3O = rate of respiration per 4 
unit of zooplankton (day ~) 
X3 = standing stock of 9 
zooplankton caeCym )., 
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The equation for the respiration rate is: 


RHO30 = RHOZ3 * exp(R*TEMP) (30) 
RHO30 = rate of respiration per | 1 
unit of zooplankton (day) 
RHOZ3 = rate of respiration at O°C. (day™*) 
R = rate of change in 
respiration rate as a 1 
function of temperature Ce. 
TEMP = mixed layer temperature ou Oy ms 


The loss of biomass by excretion is: 


VOID = TAU30 * X3 (36) 
VOID = flux of voided biomass (g C/m?/day) 


TAU30 


rate of excretion per 
unit of zooplankton -1 
biomass (day ~~) 


X3 = standing stock of 5 
zooplankton Comey im.) 


Based on data collected by Harvey et. al. (1935) 
and his own observations, Riley (1947) concluded that 
utilization (ratio of assimilation to ingestion) of grazed 
phytoplankton was very nearly unity when algal concentrations 
were low. During such periods, therefore, minimal carbon 
biomass would be voided. However, with the onset of the 
spring flowering, more food was available than could be 
used, so that quantities were voided in a semidigested 
condition. To approximate these observations, the fraction 
of biomass voided has been treated as a function of the 
grazing rate: 


TAU3SO = (0.8*(TAU23/MAXG)) * TAU23 (37) 
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TAU3O = rate of excretion per unit = 
of zooplankton biomass (day ) 

TAU23 = zooplankton grazing rate (day71) 

MAXG = maximum specific grazing rate (day™*). 


= 


When grazing rates are low (.025 day ~), the bulk of 
the ingested phytoplankton is assimilated (92%) and only 

8% is voided. During the bloom conditions observed in 

March 70% of the grazed biomass was assimilated with 30% 
being voided. These values are in general agreement with 
Conover's (p. 113, in Parsons and Takahashi, 1973) conclusion 
that, in general, zooplankton assimilation efficiencies 

range from 60 to 95%. 


The final term in the zooplankton equation represents 


losses of zooplankton biomass by predation and natural 


mortality. 
LOSS = LAM30 * X3 (38) 

LOSS = flux of biomass due to 9 
predation and mortality (g C/m” /day) 

LAM3O = rate of loss per unit of -1 
zooplankton (day ~) 

x = standing stock of 9 
zooplankton (g C/m”), 


The rate of loss has been estimated by this author as a 


constant 1% of the standing stock per day. 
LAM30 = L (39) 


L = fraction of zooplankton lost “1 
per day (day) 
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NAME 


DX1 


EKL 


FAVL 


KG 


RHOZ2 


RHOZ3 


XIMIN 


X2MIN 


TABLE I 


VALUES FOR EQUATION CONSTANTS 


AND THEIR RESPECTIVE SOURCES 


DEFINITION 


MAXN at O°C 
phosphate:carbon ratio 


nutrient concentration 


depth of Ekman layer 
nutrient regeneration 
coefficient 

half saturation constant 
for zooplankton grazing 
half saturation constant 
for phytoplankton grazing 
zooplankton loss rate 
coefficient of exchange 
maximum grazing rate 
photosynthetic constant 
respiration coefficient 


phytoplankton respiration 
veo ©. 


zooplankton respiration 
ages OoC.. 


rate of change of MAXN 
with temperature 


minimum nutrient uptake 
concentration 


minimum phytoplankton 
Concenrration for grazing 
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VALUE 





0,50 
0,774 


3,0 


00,0 


0,8 


80.0 


0.01 
0.12 
0.25 
2.9 


0.069 


0,0175 


0.019 


ails: 


G0 


Rie bene 


estimate 
Riley (1965) 


estimated from 
data 


Smith (1968) 


Butler, et al,, 
(1970) 


McAllister (1970) 


Ichimura, (in 
Parsons and 
Takahashi, 1973) 


estimate 

estimate 
McAllister (1970) 
Riley (1946) 


Riley (1946) 


Riley (1946) 


Riley (1947) 


estimate 


Dugdale (1967) 


McAllister (1970) 





IV, PROCEDURE 


A. SIMULATION 

The system of differential equations was integrated 
numerically using the System/360 Continuous Simulation Model- 
ing Program (S/360 CSMP) available on the Naval Postgraduate 
School IBM-360 computer. A brief discussion of S/360 CSMP 
is given in the Appendix, A more complete description is 


available in the IBM User's Manual (1970). 


B. DATA COLLECTION 

Theoretical results, generated from observed environmental 
fememetions (i,e€., radiation, temperature, etc,)}, were 
compared with observed seasonal variations in phosphate 
concentration and zooplankton biomass carbon. Data was 
obtained on a series of cruises conducted aboard the R/V 
ACANIA from February through November, 1974, Dissolved 
reactive phosphate was measured uSing a Technicon II Auto- 
analyser. Concentrations in the mixed layer were determined 
by computing a "weighted mean" of measured values within the 
layer. The "weight" for each value was assigned according to 
the interval between sample depths. 

Zooplankton samples were collected in a number three 
(333 um mesh) net with a mouth diameter of c.a. 25 cm. The 
carbon content of each sample was determined using a high 
temperature dry combustion thermal conductivity method 


(Traganza, Radney, and Grahm, 1975). A complete listing 
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of cruise data including phosphate, temperature, and carbon 


values is available, (Traganza, 1975). 


C, PARAMETER TESTING 

The response of the model to different environmental 
conditions was examined by generating simulations using 
hypothetical values for environmental functions. The values 
of incident radiation (RADI), mixed layer temperature (TEMP), 
mixed layer depth (Z) and upwelling velocity (W) were varied 
10% above and below their "standard value" (the observed or 
literature value). A seperate Simulation was computed for 
each variation, leading to eight seperate cases. The 
variations tested are indicated in figure (4). The O's 
represent the standard or observed value, the +'s indicate 


a 10% increase, and the X's a 10% decrease. 
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V. RESULTS 


A, COMPARISON OF SIMULATION WITH OBSERVED DATA 

The initial simulation results are shown in Fig. (5), 
Theoretical cycles appear as asterisks and the observed 
cruise values are indicated by X's. The average difference 
between Simulated and observed nutrient concentrations is 
66.4%. A 94% average error exists between the simulated 
and observed standing stock of zooplankton. Although the 
absolute error is large, the trends of the data and simulation 
are quite Similar. The principle source of error in the 
phosphate cycle was a 20 day variation between the simulated 
and observed nutrient minimum. The zooplankton simulation 
error was caused by the failure of the model to match the 
apparently low biomass values observed in late March and 
early September. 

Seasonal changes in mean phosphate concentration in the 
mixed layer were characterized by a steady increase during 
winter and spring months from an initial value of 1.5 wug-atom/1 
to a maximum 1.86 yg-atom/1. The simulated concentration 
was aSSigned an initial condition of 1.5 yg-atom/1 and subse- 
Gently increased to 1.99 uwe-atom/1. In June a sharp decline 
in concentration occurred both in the data and in the simula- 
tion. A minimum value of 0.33 yg-atom/1 was observed in early 
August. The simulated minimum (0.61 yg-atom/1) occurred in 


late August. Both the theoretical and the measured 
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concentrations increased steadily from August through 
December culminating at 1,0 ug-atom/1. 

The differences between observed and Simulated maximum 
and minimum nutrient concentrations indicated that the 
function controlling nutrient uptake, NTLIM, (eq. 27) did 
not fit existing conditions. Testing various values of the 
constants KN, (eq. 27) S and B (eq. 28) indicated that a 
reduction in the value of KN from 0,6 ug-atom/1 to 0,3 
ug-atom/1 produced a more accurate correspondence between 
measured and theoretical maximum and minimum phosphate 
concentrations (Fig. 6). 

Phytoplankton data was not available for comparison with 
theoretical values. However, the sharp decline in nutrient 
concentrations during June and July Suggests the presence 
of substantial phytoplankton stocks during this period. The 
Simulated bloom, which peaks in July, correlates with the 
observed decline in phosphate. The phytoplankton maximum 
does not, however, precede the initial zooplankton peak as 
is normally observed. 

The cycle of zooplankton biomass was characterized by 
two distinct maxima. The first peak (1.05 eC /m-) was 
observed in late July. A Simulated maximum (1.02 gC /m) 
occurred during the same period. This initial peak was 
followed by a Sharp decline of measured biomass. A minimum 
on 0.18 eC /m@ was observed in early September. A lesser 
reduction occurred in the theoretical values (0.78 eC/m-). A 


rapid increase in the standing stock followed the September 
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minimum. An overall maximum of 1,895 gC /m” was observed in 
mid-December. Simulation values also increased producing 
a lesser maximum of 1,58 gC/m?. 

The simulated energetic balance observed in the zooplankton 
varied with changes in phytoplankton concentration and 
temperature, During the onset of the phytoplankton bloom 
in mid-March, zooplankton assimilated 70% of the ingested 
carbon biomass and voided 30%, Of that carbon which was 
assimilated 60% was lost through respiration. In mid-July, 
when phytoplankton stocks had peaked, the high algal 
concentration resulted in near maximum grazing rates. Because 
more biomass was being ingested than could be used, 77% was 
voided and 23% assimilated. Due to higher mixed layer 
temperatures, 75% of the 23% assimilated biomass was required 
for respiration. 

Calculation of the annual average of monthly zooplankton 
gains and losses indicated 65% of the ingested phytoplankton 
biomass was assimilated, 35% was voided. Of the 65% 
assimilated biomass, 86% was respired leaving 14% (9% of 
that assimilated) for growth. Of this remaining biomass, 


75% was lost to predation and natural mortality. 


B. EFFECTS OF VARIATIONS IN ENVIRONMENTAL PARAMETERS 
Additional Simulations were generated to determine the 

Mesponse Of state variables to changes in the environmental 

parameters. Table II lists the percent change each modified 


parameter produced in the phytoplankton and zooplankton 
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TABLE II 


CHANGE IN MAXIMUM BIOMASS OBSERVED 


WITH +10% CHANGE IN FORCING FUNCTIONS 


PHYTOPLANKTON ZOOPLANKTON 


RADI, -10% 
RADI, +10% 
LAYER, -10% 
LAYER, +10% 
UPWEL, -10% 
UPWEL, +10% 
TEMP, -10% 


TEMP, +10% 
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standing stock maxima. Complete simulation results are 
given in Figures (7 - 10). The 10% increase in incident 
radiation produced a nearly equivalent gain in the maximum 
biomass of both phytoplankton and zooplankton (Fig. 7). 
Increased radiation produced higher rates of photosynthesis 
and larger stocks of phytoplankton. These phytoplankton 
stocks in turn supported increased grazing by zooplankton. 

Reducing the layer depth produced a slight increase in 
biomass. Confining the plants in shallower water apparently 
exposed the phytoplankton stock to higher average radiation 
intensities and increased algal availability to zooplankton. 

Surprisingly, a 10% change in upwelling velocity had 
negligible effects on plankton biomass. Increasing velocities 
did produce higher nutrient concentrations (Fig. 9). However, 
this change in concentration did not yield a Similar increase 
in production. This result stems from the non-linear relation- 
ship between nutrient concentration and photosynthetic rate 
(Fig. 2). Above 1.0 ug-atom/1, changes in nutrient concen- 
tration have a small effect on photosynthesis. Since the 
variations in the simulated nutrient concentrations were 
observed principally at concentrations above 1.0 ug-atom/1, 
these fluctuations had a minimal effect on plankton production, 
The model's response here agrees with Odum's (1971) contention 
miaee nutrients are not a limiting factor for phytoplankton 
production in regions of upwelling, 

The largest variation in zooplankton biomass (-65 to +85% 


fee the result of a +10% variation in the mixed layer 





temperature. Temperature changes directly affect the 
respiration rate of the zooplankton. Increased temperatures 
produced higher respiration rates and a net loss of biomass. 
Conversely, lower temperatures produced lower respiration 
rates and greater standing stocks. Although the change in 
daily respiration losses is only Slightly affected by temper- 
ature (7%/°C @ 15°C), the cumulative affect over a period of 


months appears Significant (see Fig, 10), 





VI, DISCUSSION 


The primary environmental factor influencing simulated 
phytoplankton stocks was incident radiation, followed 
closely by temperature effects (see Table II). The greatest 
fluctuations in zooplankton biomass resulted from changes in 
the temperature regime. The response of zooplankton stocks 
to different thermal conditions suggests a hypothesis 
concerning the rapid September (warm water period) decrease 
in zooplankton stocks, 

An examination of mixed layer temperatures and predicted 
upwelling velocities for July and August (Fig, 4) indicate 
that this period was characterized by a reduction in upwelling 
(less cold deep water moving towards the surface) which 
resulted in higher mixed layer temperatures. These higher 
temperatures may have increased zooplankton respiration 
sufficiently to produce a net decrease in biomass. In late 
fall, when lower temperatures returned with algal stocks 
remaining sufficiently large to permit grazing, a rapid 
zooplankton increase occurred. 

Certainly, other processes contribute to changes in zoo- 
plankton stocks. The late summer decline may have been the 
result of intense predation during this period. Since the 
model assumes a constant loss to predation and natural 


mortality, variations in predation have not been considered. 
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A second possibility is that water masses containing 
different size plankton stocks may have appeared in Monterey 
Bay in late summer and winter, The observed changes in the 
bay have traditionally been described in terms of three 
@ceanographic periods (Skogsberg, 1936; cited in Sverdrup 
et al., 1942). The “Davidson Current period" (relatively 
warm water) is normally observed in winter, the "upwelling 
period" (cold water phase) from early spring through summer, 
and the "oceanic period" (warm water phase), from late summer 
through fall, A comparison of these dates with the observed 
zooplankton fluctuations indicates that the first maximum 
might be associated with the upwelling period, the minimum 
with the oceanic, and the second maximum with the Davidson 
Current period. 

Although changing current patterns clearly exist in 
Monterey Bay, the role of advection in controlling plankton 
populations may be negiected if (for the time and space scale 
of interest) the region can be assumed horizontally homogeneous. 
The total rate of change of a variable "S'" as given by the 


material derivative (Se) 1s! 


one rade ae Be TPN Sores (1) 
The total material derivative (=>) will equal a simple time 
derivative ($3) DteomiorsZzonval velocities (u, v) are zero 
or if the region is horizontally homogeneous, i.e., the 


Spatial derivatives (S5, a=) are zero. 
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The argument has been made by Bolin and Abbott (1962) 
that the hydrographic and biological features observed in 
Menterey Bay are characteristic of a much wider geographic 
region. These investigators cited evidence that upwelling 
in the bay is paralleled by upwelling along most of the West 
Coast of the United States north of Point Conception. Their 
temperature observations also corroborated closely with 
the more extensive data reported in CalCOFI Report, VII, 1960. 
In addition similar phytoplankton observations were made at 
eestaction as far distant as La Jolla. Their conclusion was 
that the assumption of horizontal homogeneity iS appropriate 


for this region, 
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VII, CONCLUSIONS 


The ecosystem model Simulated with reasonable accuracy 
the observed seasonal changes in phosphate concentration 
and zooplankton biomass. The model inputs used to generate 
the theoretical results consisted of parameters defining those 
characteristics of the physical environment which effect bio- 
meereal productivity, e.g., solar radiation, temperature, 
layer depth, and nutrient input. 

Analysis of the model's sensitivity to changes in these 
environmental parameters identified the possible importance 
of thermal conditions in regulating the growth of zooplankton 
mepmlations, These theoretical results suggest a hypothesis 
heteated tO upwelling. The importance of upwelling as a 
source of nutrients has long been acknowledged. Simulation 
results indicate that upwelling may also enhance the growth 
of zooplankton stocks by bringing colder deep water to the 
surface thereby reducing zooplankton respiration requirements, 

Of course, the model can only suggest the possible signif- 
Meenee Ol this cooling effect. Actual verification of the 
hypothesis will require a precise knowledge of zooplankton 
metabolism, as well as additional measurements of phytoplankton 
and zooplankton biomass. A more detailed simulation study 
of the effects of upwelling might best be made using a two- 
dimensional model to analyze variables as a function of 


depth as well as time, 
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nutrient recycled 


(REGEN) 


vwolded biomass and 
unrecycled nutrient 


(VOID) 


nutrient inpuf 


ox = NUT + REGEN — UPTAK 
axe = PROD — RESP2 — GRAZ 
ae = GRAZ — RESP3 — VOID— LOSS 


FIGURE 1. Ecosystem Compartment Model and Flux Equations. 
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APPENDIX A 


A BRIEF DESCRIPTION OF THE S/360 CONTIUOUS SIMULATION 


MODELING PROGRAM (CSMP) 


The S/360 CSMP program is devided into three segments, 
meentitied by the statement labels INITIAL, DYNAMIC, and 
TERMINAL. The INITIAL segment is intended for the specifi- 
cation of initial conditions and for defining constants and 
parameters. The DYNAMIC segment includes the complete 
description of the system, together with any other computa- 
tions desired during the run. The TERMINAL segment is used 
for those computations desired after the completion of each 
run. The INITIAL and TERMINAL segments are optional, but 


the DYNAMIC segment is mandatory. 


structure Statements 

ptructure Statements define the model to be simulated. 
In general, rules for structure statements follow those used 
in FORTRAN. A library of special CSMP functions is also 
available. One such function used in the present model is 
mem iNTGRL function. This function computes numberically the 
integral of a function "X,"' It is expressed in the general 


form: 


¥ = INTGRL (IC, Xx) 


where the initial condition is: 
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EG 2=<¥ CO), 
This CSMP statement is equivalent to; 


t 


Y= | Sees te eC 


O 


Data Statements 

Data statements are used to assign numerical values 
to constants, parameters, and initial conditions. These 
data cards are specified by the labels CONSTANT, PARAMETER, 
and INCON respectively. If a specific parameter is to be 
tested over a range of values a sequence of simulation runs 
may be designated by enclosing several values of the variable 


in parentheses. For example: 


PARAMETER X = (5.0, 5.1, 5.2), 
specifies three simulation runs. The value of XK will be 5.0 
in the first run, 5.1 in the second, and 5.2 in the third. 

A linear function connecting data points may be 
specified by the function generator AFGEN. The set of data 


points is identified using the FUNCTION label: 


BUNeCTIONVUPWELT = CO., ,10), (10., .50). 
The points are given in sequence, X-coordinate (time) first 
followed by the y-coordinate. The program connects the data 
points with a linear function when the AFGEN function is 


called: 


W = AFGEN(UPWELL, TIME). 
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Control.Statements 
Certain operations related to translation, execution, 
and output segments of the program are specified by special 


control statements; 


etc These three labels identify the major segments 

TERMINAL of the CSMP program. 

END This statement marks the completion of the model's 
structural description, 

SORT These labels determine whether a sequence of cards 

NOSORT is to be machine sorted into a correct sequence 
or executed in the order given. 

STOP This card follows the last END statement in the 
program. 

ENDJOB Titemctauemecmuaenotes the end of a®job and must 
follow the stop card. 

TIMER This label is used with the following CSMP 
specified variables. 

PRDEL Output print increment. 

OUTDEL THine—Olot Olleput Print increment. 

FINTIM The maximum time for Simulation. 


This is only a partial listing of the functions and 
labels used in the 8/360 CSMP language. For a complete 


description see the IBM User's Manual (1970). 


67 





(AWIL SYUVIOSINSDSV=edWODN 
(GUE 8 = SSE) at) Gye) 
tae Te fs Wo Le ) iC te to oe) “(SCC 2 42852 ysl O82 eae) 
eee eee Ct) f(Gec" **991) “(occ 4nSt be eee esol) 
foe ec 7) ML eGT* * "Oy )) * (CC hee2 00) =e Ios NO lens 
SOCIO IOI IOI IOI Ga Gia kkk kok oki tak kok gk ok kok sk kak kk diake skats ate kg afc ak ak ote ate akc ak ok ofc akc abe ake ak ote oe ate ake ok le ak ot 
sx . °NS9DSV SYOLVUSNSD xx 
ee NOITLONAS YVSNIT 2FL AG GSNIASG 34V SNOILCNFA IWILINSWNOYIANS Ba 
XK aK HK aK ok ake ok aK aK oc a Xoo akc afc ake akc fc ak a alk ok ake akc ae ak ofc ok ake ak ake ok ac ok ok ofr ok aie fe ok ake ok ake eof ak ake ic ake ke a ok ok ak ake ole akg ok ak ole sk ocak ok ok ak 3c a 
JIWVNAQ 
(28T° *OLT°® *EST* J=ZDVIM Y2LIWVAVG 
C21 =9VSer Oot -ous ll *O ol =o Ve et Cova 
"8° FO°OT=NIWZX SO°O=NIWIX ‘SETSCG=S S6T0°0=E70HY *SLTIC°0=270HU 
mew og): O=4 O° 2=d "S22 °O=9XVKR 8 CT Oh *10°OSs) *S20=N\) 
eee 6" O= 1AVa ©O°OS= Ila *O°C=T XG *917-0=9 *S°O=0 ENVISNUD 
bo O=Cole i 0=coO1 eS ala Noa 
VWILINT 
eH BCR SK AS 2K 2 2 ARCA AK ACI HI COR OK NC CK aK aE RAR IR eS aK Oe AE 2K ROR RO eR ORK ROK AR IK ERC ea RC CO A oe 2 aK aK 
ibe *SYgLIWwVuVd CGNVY *SINYVLSNOD *“SNOTLIGNOD xx 
a WILINI] YOSd SSNIVA WIIYSWOW SNIVINOD INAWOSS TVILINI SHI xx 
oe ok ae oe oc ok ea ok 3 kc ake ie ke ake eae a ok ak aK ak a ake akc at ok ok 2 ake ok ole ake ak fe ke yc ok at ok ok ok oko oie oe ae ote ok ok ok ak ok ak ok fe oe ok ake ok aie ke ete ok see ok oe oe ok 
We ke OK aie a 2k ak ak ok ok sek ee eke ok ok ok ak oe a akc ae ok ok ae ak ak ok ok ok eas ake te ake oe ote ake ob kook oc ek akc ok ate ott ote ob ate aie aks ke ok ale akc ak le ok afe ate ak ok ak afe oe ak fe 
RK ok se ak a akc ake oe ake a Se ake a a ote ok oe ake ae te kc ok a io aK oie ake fe afc ake ak ak ok ok ok ca ak oc ae 2K ok ake otc ok oc ae aig ake oe ok 2 ae oe ok eke ae oi ake ale ake ak ok 
mk ok °(Z) HKLd3G YSAVT GSAXIW GNV ‘(dbisl) x 
ae AYNLVesdWw3al YSAVT G3SKIW *(M) ALIDOIVZA SONI W1SMdf *(IQvVH) NGILV + 
xe ~IOVY ILNSGISONI SONIIN] StslL3WVevd SHI "ONT T1V3MdN JO NOIONSY V NI x*x% 
ae ste ALTAILOAGOYNd IWOILSGTOIG SONSNIANI LVHL SYSlanVvd TIWINSKENOUIANS xx 
x ok YNACS NI SNOTLVIYVA IWANNY NO INGCNS3d30 3UV SLIASSY NOILVIANWIS +x 
x Xe SH ft "(€X) SSVWUIG NOLINV WOCZ7 SNOYDAIEYSH CNV *€2X) SSVWOIG xx 
xENOLYNVTdOLAHG *( TX) NOTIVYLINSJINGD SLVHdSOHd YSAVI A3XIW -S3ZIGVIYNVA «x 
ak JL1ViS 33SydHL NI SNOTLVALINIS IWNOSVSS SHL SSLVIAWIS W30COW STHL *x* 
Xe 2c ote ake ak 2 afc akc aI aK RE KC ARC Cafe aK aR He oe OC aK ae 2c a KK aK ae ake ae aK ake oe ake ake ke ake A ak aK ok te ake ak ate ate akc ak ate ate ote ak ake ake cake ak ake ak ae ake ate ak ake ok ak ak ac ok 
BOI IDIO ICI OI oI i ici dak aioli ak ak ak aio akak a aia goa dake keke ae a ak ok ag a seco ak ae ok 8c ok ah ae fe a ake a ea sk eae sake a tc oe aaa oe a ak 
WS3LSASO033 SONI TISMdN ASYSLNOW SULIL 
*x CG NISAS* dWS3// 
O9EdWS3 93X33 // 


68 





(U2epie= dX S=O> alee) 7 lay e=dy a 

(TX-TXO)x(2/W)=XIW 

(TX- TX) s(TH3/M = 1SMdN 

(2/2X)*°OOOT=9WZX 

999999 °0**19%*4S0°0+4+ 1748800 °0+40°0=% 

(Z/2X) xX EE E*EC=aTI 

((NIWTX-TXD4NH) /SONIWTX-TX)&NXUWEWIT TIN 

€4+dWsl *$ =NX VW 
ee CALVINDIVD JUV SAGILVNGS SLVY NI CaINCst SITGVITWA @LYaN #x 
SR aK KC KE akc aK ai 2 2 a 2K ae ke 2 aK a aK a IC a ae a aK ak IC a a akc a ae aie aI akc ae ak ak a ae aa ake akc afk ak ate ai ake ae ait ok akc akc akc ke ok ak ak ak ake 2 ae akc ake ie ak ak 

dWOD7*IDVAZ=7 

dhO0L* IV dSL=edWdal 

dviOIM*IVIM=M 


dhOJedx*x BVAD=IGVY 
BK AE AK IE AK ae AC i OK a “eK ak ake 3K te XS a Oe ak Ns BE He Ne ie Ae i a ais Hk ae Bie IK aS aie ak Ae aie aie aK ais aie ate ak ais 3k Me oie Xe oe oe OK ais alk als ais Xo aK OK ae SEK OK XE NK OE OK OK AK 
xe ok "9135 §9VA2 «x 
a 2 ‘ovdMm 3°] *HOL9IVA LNVISNOD Vo AG SNIWA G3LNdWOD B3HL INIAIdILINN ¥% 
Hee AG G3SLVYSNS9N 3ev SYSLIWVYVd IWINSWNCYIANS Ni SNG@ILV 1] avy ex 
FOIE I I IIa 2c i a aI a aK oi af ak ai ok a ake ake ake akc ak a aK ake ak ok ak aka ak ak ok ac oe ak kc at oka ol akc a ok fe akc ok ofc ag akc ok aie a ake ak a ak a 

(SWILYaAVIINS9S v=dW007 

On Cee 7 eo) ( Om a Vet Cea) 

preterteece § “bOc} (2 eh 6ST) {2 OO T°? bod be02 * 0) ey aAV 1 IN Gai ae 

(SWIL Sdh3LSINIODAVEGWOIL 

(Oe cle SOE) (1°CT$ “epee ss ace Gus Cee C7 sk 8 8G | eeelic) 

ee EPS elt Od pAb Tt OG yee ol- to (Oe Oy —dnas NG lente 

(SWI LS 3ISMdNIN39d4 V=dWODM 

(0268-6900. 2 602 Se ye.) 

PPPS EEE EOE TC) tale] ME RO Ze *.0 1.2.22! °.0 COC eee eae) 

eee §(99°C** C61). *(86°R! * COT) See oe Ged) *toc™ 2s |) 

coe bE gyetecy) S¢7G%b *Gh) *(°08*GT) (90% *O)=717EMdN NOTLONAS 


69 





EXxOEWVI=SSO1 
EXXOENVI=GIOA 
EXxXOEQHH=HEdS3AY 
EX*XECNVL=EZVYI 
CXxNZOHY=CdS3y¥ 
CXxCTAVL=GONd 
(Z/ 2X) S2TNVLXI=RENVIdN 
EX*xLENVL=ENIA9Sz 
XIW+ TAMGN=ELON 
ES SSNOLLYNOS S1VY WO8d GSINdWO) SHV S3XNTd 3H 
2 af ah a ake ae ae ak Re ake oe ate 2 a ake fe a 4 oi a ae ok ak ate af fe ate akc a ake ae ae a ate a ake ake 2 ake ale ake akc ak ak ok oe oe ake ak ofc of afc ak ok of ak afc ai fe ae ae ae ak ak a ak ak 2c 
T=OEWV I 
(dha ls ayako ez e=Ceur ad 
(dWI31L*Y )dX353%2Z70HY=O0COHY 
OENVIx (2/9) IAVA=TCAHVI 
ECNVIXCCOXVW/JECNVL )*8°O J=OCNVE 
O°O=ECNVI tororndeznviiai 
LYOSON 
CONIWCX-ONC?CXDFIYH) SONI WCX-OhWCXK XOX IW =ECNVEL 
WITIN+OVUxd=2INnvi 
OIC RIGOR IOI I I A IOI II a BOR RL kai ie ak ak kak sb ak ok ok ak Rak aie ak a ok ak oh ke 2 ne Ne ae he ak a aie ak oe 


He Xe ake ak OK 
a a "dsalAdWwaod 3axVv SNOT IVADS ALVS xx 
AGIOR II IC IOR IOI Hai oI i a II a fea ak ak ale ak ak ake ak ak ak ok aR kak 2k ok fe ate ake oi ok teak ak ak ok ag ok teak ake fe akc i ak ak ak ak ak akc ak at ake aie ok 


70 





gOrFONS 


dOILS 

QNJ 

(OIOA S£€dS3ae¥ 4Z2ZVYAONIEX LOTdLyd 
(2vo * dead. “COdd) Ge 10 lated 
(Hiveld? ONagsa “ih) DX eed ise 
O°4T=71901NO ‘f*+9E€=WILNIS YSWIL 
WoO] €X eT Os ex a 

1xYOSON 


(FG e Xe] } WON I=ex 

(100cxX%* 231) TYDINI=2x 

{LOOTX* 131) IYOSLNI=TX 
af ake ai ae 2k aie afc a a ok aie oleae ake ae ak ae afk ae a ae ale 28 ake kc at ak ae ae ac afc 2 fe ake ke ae ae ake ae ke ake fe ake oe fe eae ake ak ake a aft ak ae a oe fe okt afk ate ake ae afc 2K ate ak ak 2 ke 
=r "OaiNdwOd 3YV SSTIGVIYVA 3ALVLS SHI YO4 SANA MAN 3H} *“AVIVNIA *x 
Bf ae a ake Sk a ac aK aK 2k ak aca a 2K 2c He ak aK ak Se a 2 3 ae a a ak fe akc ae he ae ae ae se te afc eae 3 a fe he ke ak ae ae fe 2 ae ee ea fe aie ae a ate ake fe ake se oe aie ak Ok 

SSOT-CGIOA-€dS 3U-Z2VYN=1LODEX 

2VY9—- CdS S3u-dO0¥dd=1002X 

WIdDN-NAD IY FLAN=LOOTX 
ae Se ste ake ae akc 2k 2 ac ais oe cae fe a ake fe 2 fe ae ae fe aK fe ok fe 2c Hohe ks oe oe oe ak ae a ae ae 3k 2 ate She a he ee ak oe 2g oe ak ae aes ae oe hea 2 ok ate 3 2k 3k 


ee TSSAXNIS JO WNS JHL SV GALNdWOD B2YVESAIEGVIYVA SLVIS 32Ht Ni BSINVHOD xx 
HERE IE BE IK AIS A RE aK a CK aK aK aC aR a EK NE aK aC a IE SC aK a oa fe Kc a a a akc ake a OC aK aie ake Ree aE OS aS BRC a aE a ie ak a ai aK ae 


(an 





LIST OF REFERENCES 


Bakun, A., 1973, Coastal Upwelling Indices, West Coast of 
North America, 1946-71, NOAA Technical Report NMFS 
SSRF-671. 


Bolin, R.L., and Abbott, D,P. (1962), Studies on the Marine 
Climate and Phytoplankton of the Central Coast Area 
Of California, 1954-1960, California Cooperative 
Oceanic Fisheries Investigation Reports, v. 9, 


p. 23-45. 


mmer. F.t., Corner, £.0.5., Marshall, S.M., 1970, On the 
Nutrition and Metabolism of Zooplankton. VII. 
seasonal Survey on Nitrogen and Phosphorous 
Excretion by Calanus in the Clyde Sea Area, J. Mar. 
PiOmiss. UKs, V2 90, p. 525-560. 


Cushing, D.H., 1959, The Seasonal Variation in Oceanic 
Production as a Problem in Population Dynamics, J. 


Conomeempllor, Mer.) v. 23, p. 178-188. 


Dugdale, R.C., 1967, Nutrient Limitation in the Sea: Dynamics 
Identification and Significance, Limnol. and 
Oceanogr., v. 12, p. 685-695, 


Deedale, R.C., and Goering, J.J., 1967, Uptake of New and 
Regenerated Forms of Nitrogen in Primary Productivity, 
mee. sand Oceanoeyr., v.12, p. 196-206. 


Eppley, R.W., and Thomas, W.H., 1969, Comparison of Half- 
Saturation Constants for Growth and Nitrate Uptake 
Siecle planmicron mms Phycol.,, v. 9, .p. 375-379. 


meomine, R.H., 1939, The Control of Diatom Populations by 
Cao ee meets Explor, Mér,, ve 14, p. 210-227. 

meeoiner, A.C., 1937, Phosphate Production by Planktonic 
POU oe meme Leplor. Mer), v. 12, p. 144. 

mead. H.H.> and Braarud, T., 1935, A Quantitative Study of 
the Phytoplankton in the Bay of Fundy and the Gulf 
of Maine Including Observations of Hydrography, 
Chemtciiy samme luitbidity, J. Biol. Bd. Canada, 
Vee, Do 219-467. 


Mieris, HE,, 1959, The Nitrogen Cycle in Long Island Sound, 
Dione noneeeceanogr, Coll., v. 17, p. 31-65. 





V2 





Harvey, H.W., 1957, The Chemistry and Fertility of Sea 
Watemo —caeca., Cambridge University Press. 

Harvey, H-W., and others, 1935, Plankton Production and Its 
Conn Glew weit b1OlL, Ass. U. K., v. 20, p. 407-441. 


IBM Application Program GH20-0367-3, 1970, System/360 
Continuous System Modeling Program, User's Manual, 
IBM Corp, 


Ketchum, B.H. 1939, The Absorption of Phosphate and Nitrate 
by Illuminate Cultures of Nitzschia Closterium, Am, 
ieeesot., Vv. Zomp. 399-40". 


Kuenzler, E.J., 1965, Glucose-6-Phosphate Utilization by 
Marine Algae, J. Phycol., v. 1, p. 156-164. 


etka, A.J., 1925, Elements of Physical Biology, Williams and 
Wilkins (also issued as Elements of Mathematical 
proloey (1956), Dover. 


McAllister, C.D., 1970, Zooplankton Rations, Phytoplankton 
Mortality and the Estimation of Marine Production, 
meme hood Chains, Ed. Steele, J.H. University 
ot Calirornia Press, 


Mullin, M.M., 1963, Some Factors Affecting the Feeding of 
Marine Copepods of the Genus Calanus, Limnol. and 
Geeamear., ve SS, p. 239-250. 


wom, &.P,, 1971, Fundamentals of Ecology, 3rd. ed., W. B. 
saunders Co. 


Peamweons, IT.R. and Takahashi, M., 1973, Blological Oceanographic 
Processes, Pergammon Press. 


Patten, B.C., 1968, Mathematical Models of Plankton Productions, 
Int. Revue ges. Hydrobiol., v. 53, p. 357-408. 


Patten, B.C., 1971, Systems Analysis and Simulation in Ecology, 
v. 1, Academic Press. 


meoler, F.H., 1971, Feeding Rates. Zooplankton, In a Manual 
on Methods for the Assessment of Secondary Produc- 
tivity in Fresh Waters, Eds. Edmondson, W.T., and 
ienpere Gna. Blackwell Scientific, p. 2238-256. 


fee y, G.A., 1946, Factors Controlling Phytoplankton 
Populapmens on Georges Bank, j.°Mar. Res., v. 6, 
p. 94-73. 


Riley, G.A., 1947, A Theoretical Analysis of the Zooplankton 


Pope olasOo! soeorges Bank, J. "Mar. Res., V. 6, 
p. 104-113. 


fs 





Riley, G.A., 1951, Oxygen, Phosphate, and Nitrate in the 
Atlantic Ocean, Bull. Bingham Oceanogr. Coll., 
Veto: ante. 1% 


Riley, G.A., 1956, Oceanography of Long Island Sound, 1952- 
1954. II. Physical Oceanography, Bull. Bingham 
Oceanogr. Coll., v. 15, p. 15-46. 


Riley, G.A., 1965, A Mathematical Model of Regional 
Variations in Plankton, Limnol. and Oceanogr., v. 10 
(suppl.), p. R202-R215. 


Robinson, M.K., and Bauer, R.A., 1971, Atlas of Monthly Mean 
Sea Surface Temperature and Depth of the Top of the 
Thermocline-—North Pacific Ocean, Fleet Numerical 
Weather Central. 


Pyener, J.H., 1963, Georgraphic Variations in Productivity, 
in the Sea, v. 2, Ed. Hill, M.N. Interscience 
Publisher. 


Smith, R.L., 1968, Upwelling, Oceanogr. Mar. Biol. Rev., 
v. 6, p. 11-46. 


eeeele, J.H., 1958, Plant Production in the Northern North 
sea, Marine Res., Scot. Home Dept., v. 7, p. 1-36. 


Steele, J.H., 1959, The Quantitative Ecology of Marine Phyto- 
iomiaconeatol., Rey., v. 34, p. 129-158. 


Steele, J.H., 1974, The Structure of Marine Ecosystems, 
Harvard University Press. 


Steeman-Nielson, E., 1958, The Balance Between Phytoplankton 
and Zooplankton in the Sea, J. Cons. Explor. Mer., 
Veeco wep. Lf o—186. 


ewerdrup, H.U., 1953, On Conditions for the Vernal Blooming 
Steluvtoulankwvon, J. Cons. Explor. Mer., v. 18, 
p. 287-295. 


menarup, H.U., Johnson, M.W., and Fleming, R.H., 1942, 
PneCwO@etis me nelmeehycics, Chemistry, and General 
Biology, Prentice Hall, 


[mneek.v., and DeSanto, R.S., 1972, Elements of Marine 
Pero. 2a ceca, . opringer-Verlag. 


ieee anza, &.D., 1975, Investigation of Biochemical Relation- 
ships for Determining Concentrations of Zooplankton 
Biomass and its Correlation with Chemical and 
Acoustical Properties of the Ocean, ONR, code 480 
project, (unpublished data). 


74 





iain et Pee nadieyermo.C,, and Grahm, K.J., 1975, A 
Convenient High Temperature Dry Combustion Thermal 
Comuuectiviuy  nstrument for the Determination of 
Carbon in Marine Zooplankton, Naval Postgraduate 
School Technical Report No. NPS-o8Tg75081 . 


Volterra, V., 1926, Variations and Fluctuations of the Number 
of Individuals in an Animal Species Living Together, 
leeecom wrx Dlom.  Mer.,, v. 3, p. l-ol. 


Walsh, J.J., and Bass, P.B., 1971, Oceans, A Seagoing 
Simulation Program-A User's Guide to the University 
of Washington's IBM 1130 Spatial Version of COMSYS 1, 
Special Report No. 48, Dept. of Oceanography, 
University of Washington. 


Weoeter, W.s., and Reid, J.L., 1963, Eastern Boundary 


CGuereivc. im ene cea, V. 2, Ed, Hill, M.N., 
Interscience Publishers. 


aS 








LO}. 


1 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Dbibprary (Code 0212) 
Naval Postgraduate School 
Monterey, California 93940 


Department of Oceanography, Code 58 
Naval Postgraduate School 
Monterey, California 93940 


Oceanographer of the Navy 
feriman Building No. 2 

200 Stovall Street 
Alexandria, Virginia 22332 


Office of Naval Research 
Code 480 
Weernecuon, Virginia 22217 


Dr. Robert E. Stevenson 

Dementiric Liaison Office, ONR 
Scripps Institution of Oceanography 
MaeJotla, California 92037 


beprary, Code 3330 
Naval Oceanographic Office 
masmineton, D.C. 20373 


sO Library 

University of California, San Diego 
Pow Dox 2367 

foeoolla, California 92037 


Department of Oceanography Library 
University of Washington 
seattle, Washington 98105 


Department of Oceanography Library 
Oregon State University 
Corvallis, Oregon 97331 


Commanding Officer 


Fleet Numerical Weather Central 
Monterey, California 93940 


76 


No. 


Copies 





12. 


13. 


14. 


15. 


nS. 


17. 


He . 


To). 


20. 


ZA. 


22. 


Department of the Navy 

Commander Oceanographic System Pacific 
Bex 13590 

FPO San Francisco 96610 


Commanding Officer 
Pi ronment tabe:rearcrron Research Facility 
Monterey, California 93940 


CAPT Woodrow W. Reynolds, Code 58 
Department of Oceanography 

Naval Postgraduat School 
Monterey, California 93940 


LCDR Calvin Dunlap 
COMOCEANSYSPAC 
FPO San Francisco, 96601 


Mr. Andrew Bakun 

National Marine Fisheries Service 
Naval Auxiliary Landing Field 
Monterey, California 93940 


Dr. Alan R. Washburn 

Department of Operations Research 
Naval Postgraduate School 
Monterey, California 93940 


Dr. Eugene D. Traganza, Code 58 Tg 
Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93940 


Dr, Stuart E. Knapp 

Department of Veterinary Medicine 
Oregon State University 
Corvallis, Oregon 97330 


Dr. Stephens P. Tucker, Code 58 
Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 


LT Robert T. Pearson 
SWOSCOLCOM 

eo. Naval Station 
Newport, Rhode Island 


Dr. D. P. Abbott 
Hopkins Marine Station 
Poactirer Grove, California 93950 


U7 

















fa ae 


_ 


Thesis Lea 
| P3172 Pearson 
Mc. A computer simulation 


model of seasonal 
variations in ocean 
production for a region 
of upwelling. 





, _ 
TT 





: 





> 





